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A multi-objective performance optimization method is proposed, and the problem that single structural parame- 


ters of small fan balance the optimization between the static characteristics and the aerodynamic noise is solved. 


In this method, three structural parameters are selected as the optimization variables. Besides, the static pressure 


efficiency and the aerodynamic noise of the fan are regarded as the multi-objective performance. Furthermore, the 


response surface method and the entropy method are used to establish the optimization function between the op- 


timization variables and the multi-objective performances. Finally, the optimized model is found when the opti- 


mization function reaches its maximum value. Experimental data shows that the optimized model not only en- 


hances the static characteristics of the fan but also obviously reduces the noise. The results of the study will pro- 


vide some reference for the optimization of multi-objective performance of other types of rotating machinery. 
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Introduction 


Small axial flow fans have the advantages of simple 
structure and convenient installation, so they are com- 
monly used as preferred heat dissipating components in 
cooling system of electronic products, such as high pow- 
er servers etc. The performance and noise conditions of 
the fan will directly influence the working life and the 
comfort of the operation environment of electronic prod- 
ucts. Many scholars have found that it is difficult to im- 
prove the static characteristics of the fan and reduce the 
aerodynamic noise by optimizing the single structural 
parameters at the same time. Literature!!! analyzed the 
influence of the hub on the performance of the small axi- 
al fan based on the experimental and theoretical research, 
and found that the static pressure efficiency was reduced 
when the hub was too large or too small, and then pro- 
posed the optimal hub design for small axial fans with 
the aerodynamic noise of the fan as the optimization tar- 
get. Literature’*! studied the influence of the struts and 
the circumferential bending angles of the blade on the 


aerodynamic performance of the small axial fan, and 
found that the improvement of the aerodynamic perfor- 
mance of the fan was limited only by optimizing the 
bending angle of the blade or the struts. Literature!" stu- 
died the effect of tip end-plate on the noise performance 
of the fan, and it was concluded that the method of add- 
ing tip end-plate to impeller blades had a positive influ- 
ence on reducing noise, but it may diminish the static 
characteristics of small axial fan to some extent. Litera- 
turel! applied theories of isolated blade design method 
and genetic algorithm to optimize the hub ratio of axial 
flow fan. The optimization objective was the flow of the 
fan, and the results showed that there were varied im- 
provements on the static pressure after optimization. Li- 
terature!®! selected the straight blade of a small axial flow 
fan as the research object, set the static pressure as op- 
timal objective and used genetic algorithm and approx- 
imate function method as the optimization methods. The 
better performance of the fan model by optimization cal- 
culation was obtained and the static pressure rose by 
nearly 20%. The single target performance optimization 


of the fan revolved in the above literature, and the algo- 
rithm for multi-objective performance optimization was 
rarely mentioned in these literatures. 

Relying on the previous research of the project group, 
a multi-objective optimization method is proposed based 
on the response surface method and the entropy method, 
and the aerodynamic performance of the optimized mod- 
el is carried out. This method selects several structural 
parameters as the optimization variables which have 
large influence on the static characteristics of the fan, and 
selects the static pressure efficiency and the noise of the 
fan as the multi-objective performance. 


Selection of optimization parameters 


There are many structural parameters that affect the per- 
formance of the fan. In this paper, the parameters such as 
blade number", stagger angle'*!, hub ratio”!, tip clearance!” 
are discussed and the influences of the above parameters on 
the static characteristics and the internal flow field of the fan 
are analyzed in the previous work. The main structural pa- 
rameters that affect the static characteristics of the fan are 
obtained from four structural parameters by orthogonal test. 
Blade number, hub ratio, tip clearance are chosen as the 
optimization variables. The original data for these optimiza- 
tion variables is taken from the fan model M1, which is 
derived from a small axial fan with heat dissipation 


Fig. 1 The appearance of the fan M1 
within a certain brand server. The working voltage of the 
fan is 12 V DC and the rated speed of the fan is 3000 
r/min. The appearance of the fan is defined in Fig.1, and fan 
structural parameters are listed in Table 1. 


Table 1 Main parameters of the fan M1 


Mai parameters Value 
External diameter(mm) 84 
Casing height(mm) 24 
Blade number 7 
Tip clearance(mm) 1.5 


Hub ratio 0.40 


Stagger angle(°) 46.9 
Blade gap(mm) 15 
Blade width(mm) Root 16/tip 30 


Multi-objective optimization method 


Response surface model 


The response surface algorithm based on Box-Behnken 
design is an optimization method which not only can use 
the visual graph to reflect the relationship between the 
response target and the factor, but also give the function 
relation by modeling. The modeling is based on the ex- 
perimental data, and the first order response surface 
model or the second order response surface model is of- 
ten used to construct the response surface. The second 
order response surface model applied in this paper is: 


m m 


y=B, t $ px, t Bix; H È Ras (1) 
i=l i=l 


i<j 


Where x; is the design factor, p; indicates the linear effect 
of x; Piu is the second order effect of x;, and J; represents 
the interaction effect of x; and of x;. Three factors includ- 
ing blade number, hub ratio and tip clearance are used in 
the response surface model, while two targets such as 
static pressure efficiency and noise are needed to be re- 
sponded. In this paper, the function relation among three 
structural parameters of the fan and the static pressure 
efficiency and noise of the fan is established by the re- 
sponse surface algorithm, and the optimal solution is 
obtained by the function relation. 

Before optimizing the response surface model, it is 
necessary to first determine the range of design factors 
according to the numerical results of the single structural 
parameter. The feasibility of the numerical calculation 
has been demonstrated in the literature!””!, which this 
article will not repeat. The range of blade number is from 
7 to 11, and the range of tip clearance is from 1mm to 
2 mm, and the range of hub ratio is from 0.30 to 0.49. 
Each factor selects three levels to meet Box-Behnken 
design requirements for second order response surfaces. 
Taking into account the normative nature of the test, 
three structural parameters above are represented with 
design factors X;, X2 and X3. The different levels of three 
design factors are coded. The levels of each factor and 
coding are shown in Table 2. 


Table 2 The levels of the factors in response surface model 


The levels and the code 


Factors 
-1 0 1 
Blade number 7 9 11 
Tip clearance/ mm 1 1.5 2 
Hub ratio 0.30 0.40 0.49 


According to the Box-Behnken design requirements of 
the response surface algorithm, the test program requires 
a total of 12 sets of the combination of different structur- 
al parameters of three-dimensional models and 5 sets of 
original models to complete the test. The data of static 
pressure efficiency and noise of 17 different fan models 
is obtained by the test. The efficiency is obtained at the 
rated flow rate of 0.011 kg/s, and the noise is got at the 
position where is at 1m away from the fan along the axial 
direction. The test protocol and test results are shown in 
Table 3. 


Table 3 The test protocol and test results 
No. X X,/mm X3 Efficiency/%  Noise/dB 
1 -1(7) 001.5) 00.40) 31.5 40.53 
2 -1(7) 1(2)  0(0.40) 28.35 38.94 
3 o0)  -1(1)  1(0.49) 29.23 42.58 
4  -(7) 001.5) -1(0.30) 27.57 39.81 
5 109 1(2)  0(0.40) 28.56 43.23 
6  -() 0(1.5)  1(0.49) 28.19 40.24 
7 -() 01.5)  0(0.40) 31.57 40.52 
8 0(9) 1(2) -1(0.49) 29.62 41.93 
9 -1(7)—«0(.5)—«0(0.40) 31.55 40.54 
9 111) 0(1.5) — -1(0.30) 28.81 42.73 
wu -7)  -1(1)  0(0.40) 27.96 39.64 
2- -I(7) 001.5)  0(0.40) 31.46 40.51 
13 o)  -1(1)  -1(0.30) 28.62 41.53 
4 11) 001.5) 1(0.49) 29.45 43.01 
15 11) -I(1)  0(0.40) 29.12 43.51 
6 0) 1(2)  -1(0.30) 28.93 41.32 
17 -1(7) 001.5)  0(0.40) 31.48 40.49 


Analysis of test results 


The multivariate linear regression model of the static 
pressure efficiency of the fan is established by the re- 
sponse value of the test data in Table 3, as shown in equ- 
ation (2). 

Y, = -43.27592 +8.7091 1X, +16.62068X, 
+107.73850X, —0.23750.X,X, +0.026316X,X, (2) 
+0.42105X,X, —0.45119X,7 —4.83900X,,” 
-133.21330 X, 
Where Yj is the response value of the static pressure effi- 
ciency of the fan, X, indicates blade number, X> is tip 
clearance and X; represents hub ratio. 

Equation (2) is calculated and found in the formula 
that the maximum value of Y; is 31.5661 when X; is 9, X% 
is 1.51 mm and X; is 0.41. Meanwhile, the minimum 
value of Y, is 26.2128 when X; is 7, X is Imm and X; is 
0.30. 

Similarly, the multivariate linear regression model of 
the aerodynamic noise level of the fan is established by 


the response value of the test data in Table 3, as shown in 
equation (3). 

Y, = 53.64014- 0.15779X, —7.71726X, —54.63463.X, 

+0.10500X,X, -0.19737 X, X, -2.31579 X-X, (3) 
+0.052437 X" + 2.40900X,° + 79.75069.X, 
Where Y, is the response value of the aerodynamic noise 
level of the fan, X, indicates blade number, X> is tip 
clearance and X; represents hub ratio. The calculation 
results show that the minimum value of Y> is 38.9178 dB 
when X; is 7, X% is 1.63mm and_X3 is 0.38. Meanwhile, 
the maximum value of Y> is 44.2737 dB when X; is 11, 
X is 1mm and X; is 0.49. 


Establishment of multi-objective optimization model 


In this paper, high efficiency and low noise are the 
multi-objective optimization objectives for fan perfor- 
mance. Since the dimension of the static pressure effi- 
ciency of the fan is different from that of fan noise and 
they are different orders of magnitude, the static pressure 
efficiency and noise of the fan are prepared by the nor- 
malization method. The normalization method can trans- 
form the different indexes into the function coefficients 
without dimension, and then establish the comprehensive 
evaluation function according to the weight coefficient of 
each index. The static pressure efficiency ¥ can be di- 
mensionless and the value corresponds to the interval 
from 0 to 1 by the equation (4). Similarly, the aerody- 
namic noise level Y, can be dimensionless by the formula 
(5) and the value corresponds to the interval from 0 to 1. 
The multi-objective optimization function F is estab- 
lished on the basis of equation (4) and equation (5), as 
shown in equation (6), where a and b denote the weight 
coefficient of the static efficiency and noise, respectively. 
The multi-objective optimization problem is transformed 
to the single-objective problem of function F by equation 
(6). The target with high efficiency and low noise can be 
transformed into the maximum value of function F. 

Y, = Yimin 


dl= (4) 
Yi max = Yun 
d2- ux- h (5) 
Y, max ~ Y, um 
F =ad1+bd2 (6) 


The weight coefficients a and b of the function F are 
determined by the entropy method which is an objective 
weighting method. The weight of the index in the entropy 
method is obtained by calculating the information entro- 
py which is correlated to the change degree of the index 
on the whole system. The basic steps of solving the 
weight coefficients a and b of function F using the en- 
tropy method are as follows: 

(1)The different indicators need to be standardized 
according to the equation (7) or equation (8). 
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i X; — X min 
Kip 
X max— X min (7) 
g X max— X; 
X al 
X max— X min (8) 


Where x; denotes the value of the i indicator, x, is the 
maximum value of the 7 indicator, Xmin is the minimum 
value of the i indicator, and x’; is the standardized value. 

(2)The normalized values are scaled to the specific 
weight P; according to the equation (9). 


p=! (0< P;<1) 


ij m 
F3 (9) 


i=l 
(3) The entropy e; of each index is determined, and the 
expression is: 


e=-K}_ P, In Pr (10) 
where K is the constant that is greater than zero, and K 
1 
generally takes ——. 
Inm 


(4)The differential coefficient d; of each index is 
solved. The differential coefficient d; denotes the differ- 
ence between the entropy e; of each index and the value 
of 1, and its value directly affects the proportion of the 
weight. The greater the differential coefficient is, the 
more important indicator it is. The expression is: 


dj=1-e, (11) 


(5)The weight value w; of each index is determined, 
and the expression is: 
dj 


Fa (12) 


According to the above steps, the dimension data of 
each index can be used to calculate the proportion of 
each index. In this paper, the index of static efficiency 
and the noise of the fan are taken from the experimen- 
tal data of the second order response surface in Table 3. 
According to the basic steps of the above entropy me- 
thod, the weight coefficient of static pressure and noise 
can be determined. The weight coefficient a is 0.5240 
and b is 0.4760. The weight coefficients a and b are 
substituted into the equation (6) to determine the ex- 
pression of the multi-objective optimization function F 
as follows: 


Wj = 


F =0.5240d1+ 0.4760d2 (13) 

Substituting equations (4) and (5) into equation (13), 

the mapping relationship between the function F and 

the structural parameters including X;, X, and X; can be 
established, as shown in equation (14). 


F = -7.6342 + 0.8665 X, + 2.3128X, +15.4014.X, 
—0.03255X X, +0.24701X,X, + 0.02X,X, (14) 
—0.04886 X°, — 0.6878,” —20.1272X,° 
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According to the equation (14), the maximum value 
Fax Of the function F is 0.8641 when X; is 9, X, is 1.56 
mm and X; is 0.39. The model constituted by the above 
structural parameters is redesigned and finally the opti- 
mized fan model M2 is formed, as shown in Fig.2. 


Fig. 2 The optimized fan model M2 
Performance analysis of optimized fan model 


Static characteristics 


Fig.3 shows the experimental curve of static characte- 
ristics of the optimized model and the prototype model. 
The experimental test process in the literature [4] has 
been analyzed in detail, so it is no longer discussed sepa- 
rately in this article. It can be seen from the figure that 
the static pressure of the fan model M2 is obviously larg- 
er than that of the prototype fan in the whole flow range. 
When the flow rate is less than 0.005 kg/s, the difference 
of the static pressure between two fans is small. When 
the flow rate gradually increases, the difference of the 
static pressure between two fans becomes larger than 
ever, besides the maximum difference is about 6 Pa. The 
efficiency of the optimized model is larger than that of 
the prototype in the whole range, besides the efficiency 
difference is 5.81% in operating condition. 
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Fig. 3 Static characteristics of the fans 


Aerodynamic noise characteristics 


The experimental curve of the power spectral density 
of aerodynamic noise for the model M2 and M1 is shown 
in Fig.4. The monitoring point is 1m away from the fan, 
which is located in the far field propagation region of the 
aerodynamic noise. It can be seen from the figure that the 
power spectral density of aerodynamic noise of the pro- 
totype model is mainly distributed in the low frequency 
band within 200Hz and the amplitude of the power spec- 
tral density is strongly pulsed. However, the power spec- 
tral density distribution of the optimized model is uni- 
form and the amplitude is very small. Obviously, the 
power spectral density of the optimized model in low 
frequency band is smaller than that of the prototype 
model. The noise at 1m away from the fan is mainly the 
broadband noise. 

Fig.5 shows the experimental curve of the noise spec- 
trum characteristic of the optimized model and the pro- 
totype model at Im away from the fan. It can be seen 
from the figure that the noise reduction effect is very 
obvious since the sound pressure level of the optimized 
model in the range of most frequency bands except 200 


Hz-300 Hz is smaller than that of the prototype model in 
the rated flow condition. The aerodynamic noise level of 
the optimized model is 5.8 dB lower than that of the pro- 
totype model. 
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Fig. 4 The power spectral density distribution of two fan 
models at 1m away from the fan 
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Fig. 5 Spectral characteristic curve of two fan models at 1m 
far away from the fan 


Internal flow field analysis 


Fig.6 shows the experimental curve of the velocity of 
the prototype model and the optimized model along the 
blade height direction near the outlet of the fan. It can be 
seen from the figure that the variation trend of the radial 
velocity of the prototype model and the optimized model 
along the blade height direction is basically consistent. 
The radial velocity distribution of the optimized model in 
the mainstream region (blade height 0-35mm position) is 
more stable than that of the prototype model, besides the 
amplitude is also smaller, which is conducive to the fan 
energy transmission and enhances the fan transmission 
efficiency. There is only one point of the positive and 
negative fluctuations in the radial velocity of the opti- 
mized model, which is only distributed in the middle of 
both positive and negative fluctuation points of the radial 
velocity of the prototype model. It indicates that the vor- 


tex flow condition in the region is improved and the 
aerodynamic noise reduces. The axial velocity of the op- 
timized model is similar to that of the radial velocity in 
the mainstream region. The maximum velocity of the fan 
is distributed in the tip region of the fan, and the maxi- 
mum axial velocity of the optimized model is 16.7% 
larger than that of the prototype model. The high speed 
flow in the tip area of the fan is beneficial to overcome 
the influence of the boundary layer, so as to improve the 
flow condition of the tip of the fan. 
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(a) The radial velocity curve 
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(b) The axial velocity curve 


Fig. 6 The velocity curve along the blade height direction 
near the outlet of the fan 


Fig.7 shows the pressure pulsation in the middle of the 
clearance when the working condition is in the unsteady 
flow condition. As can be seen from the figure, the pres- 
sure value of the optimized model is generally higher 
than that of the prototype model. The pressure pulsation 
of the optimized model is generally between 40 Pa and 
80 Pa, while the pressure pulsation of the prototype mod- 
el is mainly between -120 Pa and 5 Pa. It is clear that the 
value of the pressure pulsation of the prototype model is 
significantly higher than that of the optimized model, 
which inevitably confirms that the aerodynamic noise 


characteristics of the optimized model are better than 
those of the prototype model. 
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Fig. 7 The pressure pulsation in the middle of the clearance 


Based on the above analysis of the aerodynamic per- 
formance of the optimized model, it can be concluded 
that the performance optimization method adopted in this 
chapter is effective, which provides a kind of method and 
reference for other types of fan performance optimiza- 
tion. 


Conclusion 


In this paper, the multi-objective optimization of static 
pressure and aerodynamic noise of the fan is carried out 
based on the structural parameters regarded as the opti- 
mization variables. The experimental results show that 
the static characteristics and aerodynamic noise characte- 
ristics of the optimized model are better than those of the 
prototype model in the rated flow condition. The overall 
static pressure of the fan is increased by 6Pa, and the 
efficiency of the fan is improved by 5.81%, besides the 
aerodynamic noise of the fan in the far field reduces by 
5.8dB. The feasibility of the multi-objective optimization 
strategy used in this paper is verified by the experimental 
results. 


Acknowledgement 


This work was supported by Open Foundation of Zhe- 
jiang Provincial Top Key Academic Discipline of Me- 
chanical Engineering and Zhejiang Sci-Tech University 
Key Laboratory (ZSTUME 01A04). 


References 
[1] Grimes R, Quin D, Walsh E, et al. Theoretical and expe- 


rimental investigation of the scaling of micro fan perfor- 
mance[C].ASME International Mechanical Engineering 


[2] 


[B] 


[4] 


[6] 


Congress and Exhibition, Washington D.C, USA, 2003, 
24(2): 63-69. 
Chu Wei. Research on influence of struts on aerodynamic 


performance noise of small axial flow fans [D].Hangzhou: 


Zhejiang Sci-tech University Master’s Thesis, 2015. 
Chu Wei, Jin Yingzi, Wang Yanping. Effect of Circumfe- 
rential Skewed Blades on Static Performances and Inter- 
nal Flow Characteristics of Small Axial Fans [J].Journal 
of Zhejiang Sci-tech University(Natural Sciences), Vol. 
33(6), pp. 818-823, (2015). 

Mao Hongya, Wang Yanping, Lin Peifeng, Jin Yingzi. 
Influence of Tip End-plate on Noise of Small Axial 
Fan[J]. Journal of Thermal Science, Vol.26(1), pp.30-37, 
(2017). 

Zhao Yi, Jin Yingzi, Wu Wenhao, Zhang Li. Small Axial 
Fan Structure Optimization Research Based on the Anal- 
ysis of the Internal Flow Fan{J]. Journal of Zhejiang 
Sci-Tech University, Vol.29(1), pp. 89-93, (2012). 

Zhu Lifu, Jin Yingzi, Zhao Yi, Wu Yongmin. Study on the 
Optinization of Small Axial Fan Based on the Circumfe- 
rence of Blade Profile Parameterization[J]. Journal of 


[7] 


[8] 


[10] 


Zhejiang Sci-Tech University, 
(2013). 
Zhang Li, Jin Yingzi, Wang Yanping, Wu Chuanyu. Ex- 


Vol.30(1), pp.81-86, 


perimental and Numerical Investigation of Blades Effect 
on Aerodynamic Performance of Small Axial Flow Fan. 
Proceedings of 3rd Asian Joint Workshop on Thermo- 
physics and Fluid Science , Matsue, Japan,Sept 10-13, 
2011, 87-91. 

Lihong Wang, Yingzi Jin, Baoling Cui, Yuzhen Jin, Jin 
Lin, Yanping Wang, Chuanyu Wu.Factors affecting small 
axial cooling fan performance. Journal of Thermal 
Science, Vol.19(2), pp.126—131, (2010). 

Pin Liu, Yingzi Jin, Yanping Wang, Baoling Cui. Effect 
of hub-ratio on performance of small axial flow fan. Pro- 
ceedings of 2010 International Conference on Pumps 
and Fans, Hangzhou, China, October 18-21, 2010, 
(B006): 1-5. 

Wang Lihong, Jin Yingzi,Zhang Li. Numerical Simula- 
tion on the Small Cooling Fan Internal Flowing [J]. 
Journal of Zhejiang Sci-Tech University, Vol.28(1), 
pp.73-78, (2011). 


